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ABSTRACT 


This study offers a novel solution to deal with the low signal-to-noise ratio 
and slow execution rate of the first derivative edge detection algorithms 
namely, Roberts, Prewitt and Sobel algorithms. Since the two problems 
are brought about by the complex mathematical operations being used by 
the algorithms, these were replaced by a discriminant. The developed 
discriminant, equivalent to the product of total difference and intensity divided 
by the normalization values, is based on the “pixel pair formation” that 
produces optimal peak signal to noise ratio. Results of the study applying 
the discriminant for the edge detection of green coffee beans shows 
improvement in terms of peak signal to noise ratio (PSNR), mean square error 
(MSE), and execution time. It was determined that accuracy level varied 


Roberts according to the total difference of pixel values, intensity, and normalization 

Sobel values. Using the developed edge detection technique led to improvements 
in the PSNR of 2.091%, 1.16 %, and 2.47% over Sobel, Prewitt, and 
Roberts respectively. Meanwhile, improvement in the MSE was measured to 
be 13.06%, 7.48 %, and 15.31% over the three algorithms. Likewise, 
improvement in execution time was also achieved at values of 69.02%, 
67.40 %, and 65.46% over Sobel, Prewitt, and Roberts respectively. 
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1. INTRODUCTION 

Coffee has achieved notable status in the world market and is a major influencer on culture 
and economy of many nations [1]. In the current world ranking, coffee is ranked second to water as the most 
consumed food product [2] and the second most in-demand product next only to petroleum [3]. Despite this, 
there is not a single, unified coffee bean grading standard accepted and adopted by all coffee exporting 
countries [4]. Each nation either follows an existing standard or has created its own standards [5]. Before green 
coffee beans can be classified and graded, their features must first be extracted [6]. Many kinds of research 
have been done in green coffee beans feature extraction which includes the use of chemicals [7, 8], different 
types of spectroscopy such as fluorescence [9], near infrared (NIRS) [10, 11], Fourier transform [12], 
and Raman [13, 14], electronic tongue [15], electronic nose (16), and image processing [17, 18]. 

The main disadvantage of using chemicals for feature extraction is the disposal of chemicals 
after use which are oftentimes hazardous to the environment [19]. The main drawback of the use of different 


Journal homepage: http://journal.uad.ac.id/index.php/TELKOMNIKA 


2028 O ISSN: 1693-6930 


types of spectroscopy is it requires costly machines [20]. Electronic nose and electronic tongue also require 
complicated sensors that are not very easy to set-up and used [21-23]. Aside from these major disadvantages, 
the above methods also need to destroy the beans as they require grinding, require rigorous sample preparation, 
and has lengthy feature extraction in terms of time [24]. A notable exception is the use of near-infrared 
spectroscopy which does not need to destroy samples, requires very simple sample preparation, and has short 
span of time for feature extraction. Nevertheless, NIRS machines are expensive [25]. 

Image processing is one of the most promising techniques to extract green coffee bean features 
because it does not require complicated sensors and highly technical electronics [26, 27]. It only needs 
a camera, a lighting mechanism, a computer or a microcontroller and a good algorithm for feature extraction 
[28, 29]. A principal operation in image processing is the detection of edges [30]. Edges are the borderline 
in the image revealed by pixels exhibiting disruption in grey level to adjoining pixels [31]. To sustain 
the arrangement of an object while avoiding loss of information for another image processing is the primary 
objective of using edge detection [32]. 

The detection of the green coffee features using image processing require edge detection algorithms. 
This can be done using a variety of techniques of which first-order derivative-based algorithms such 
as Sobel, Prewitt, and Roberts are the simplest and most commonly used. However, these algorithms 
are susceptible to noise due to diverse modes of differential operations [33, 34]. At the same time, 
the computational complexity of the three algorithms has also resulted to large execution times [35-37]. 
As such, there is a need to address the susceptibility to noise of the first-order derivative-based algorithms 
as well as their slow execution times and test their applicability towards the edge detection and feature 
extraction of green coffee beans. 


2. RESEARCH METHODS 
2.1. Design considerations 

The major considerations for the study are the type of hardware and software to use to achieve 
the objectives of the study. It was decided that the suitable software to use is the Python because of the Open 
CV which has many libraries in image processing. Since the Raspberry Pi uses Python, it was the hardware of 
choice. The Raspberry Pi camera module was not used because an ordinary webcam can be used for 
the same result of the Raspberry Pi camera. The cost was also considered as the A4 tech webcam is half 
the price of a Raspberry Pi Camera. 


2.2. The developed new edge detection algorithm 

In this study, the differences between pixels surrounding the pixel of interest as shown in Figure 1 
was computed. The pixel of interest is the pixel in the middle or the pixel that is surrounded by 8 other pixels, 
namely the upper left, left, lower left, top, bottom, upper right, right and lower right. To calculate the total 
difference, (1) was used. 


Total Difference = (Top Pixel Value-Bottom Pixel Value)+ (1) 
(Upper Left Pixel Value-Lower Right Pixel Value)+ 
(Left Pixel Value-Right Pixel Value) + (Lower left Pixel Value-Upper Right Pixel value) 


The sample computation of total difference where the pixel of interest is an edge is shown 
in Figure 2. Using formula number 1, all values of the total difference greater than zero is considered as 
an edge, while non-edge will have a value of zero. The discriminant for the new edge detection algorithm is 
the total difference multiplied by intensity and divided by normalization value as shown in (2). 


l hy tes ee . l 
New Edge Detection Discriminant = Total Difference * ee 


Normalization (2) 

The pseudocode for edge detection is shown in Table 1. Using the program shown as pseudocode 
in Table 1, experiments were made to determine which possible combinations of pixel differences between 
pixel values of top and bottom, left and right, top left and top right bottom left and bottom right yielded 
the best edge detection. The best combination based on the output produced is the total difference of top pixel 
minus bottom pixel plus left pixel minus right pixel that is the reason why other differences were not included 
in the pseudocode. 
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Figure 1. Differences of pixel value was Figure 2. Sample computation of the total difference 


computed to classify pixel of interest 


Table 1. The pseudocode for the new edge detection algorithm 
Algorithm 
# program start 
Import OpenCV library 


start 
# input section 
define the normal value for x 
x value to be divided by 500 
define the new edge algorithm parameters - image and intensity 
define height I and width j of an image 
define the edge 
# processing section 
for all height i and width jJ pixels in range 
extract pixel values of 
top and bottom 
left and right 
top left and top right 
bottom. lert and bottom. rignt 
extract differences 


difference 1 = top minus bottom 
difference 2 = left minus right 
extract total difference 
total difference = difference 1 + difference 2 
total difference = normal (total difference) * intensity 
extract pixels of the image 
image pix = image [i , j] 
extract edge image 
edge image | i , j] = image pix * total difference 


# output section 
display input image 
display input image converted to gray scale 
display detected edge 


end 
#eop 


A normalizing value is a ratio, which means that for any value of total difference, it is divided 
by 500. The normalized value is then multiplied to the intensity which was given a value of 10. Without 
the normalizing value and intensity value, the edge that was produced by using only the total difference is not 
continuous. In the pseudocode shown in Table 1, the normalized value is 1/500 and x is the total difference. 
The program flowchart of the system is shown in Figure 3. 

The green coffee beans samples were arranged manually by hand to the sample feeder. The whole 
program was written in Python for the prototype. The start of the program is by setting up the slider button so 
that the user can set-up the Normalize value. Then the image is captured by the camera and edge is detected 
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using the newly developed algorithm. The edge is necessary to determine if the object is a coffee bean or not. 
If the object has an edge, therefore, it is a coffee, the contour will be determined and a bounding box will fence 
the coffee beans. The bounding box is the region of interest. 


2.3. The developed hardware for the prototype 

The user can access the prototype through the graphical user interface shown in Figure 4. 
The graphical user interface (GUI) of the prototype has 5 buttons corresponding to different functions, namely 
normalize, capture, edge, data, and save. The Normalize button is a slider button where the value can be 
adjusted from 0 to 100. The Capture button is for taking the images of the coffee beans. In Figure 4 there are 
two images captured by the prototype, the sizes of these images are 300x300 pixels and 800x600 pixels 
respectively. Edge button is for showing the edge detected by the developed algorithm. The Data button is for 
extracting the features of the coffee beans. When the edge is detected, the program then determines the contour 
inside the edge. A rectangular bounding box indicates that the edge and contour have been detected and inside 
the bounding box the morphology, color, and texture will be extracted. The save button is for saving the 
extracted features in the USB flash drive. The files can be viewed using the Libre Office Calc in the spreadsheet 


format as shown in Figure 5. 
Capture image 


Find edge 


Feature 
extraction 
Save feature inb 
Variable 


Figure 3. The program flowchart for the system 










Mo 








Figure 4. The graphical user interface Figure 5. The extracted features viewed using 
the Libre Office Calc 


The first column of the Libre Calc Office spreadsheet shown in Figure 5 is the sample number 
of the coffee bean and the rest of the columns are the 21 extracted features from that bean. The bottom view of 
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the prototype is shown in Figure 6. The LED strip powered by 12 V DC adapter serves as the lighting 
mechanism of the prototype shown in Figure 6. Because of the major advantage of the developed algorithm 
wherein, the image parameters can be adjusted using the Normalized button, the device can function even 
without the LED strip ON as long as there is ample light in the surrounding area or in the room. In the center 
of Figure 6 is the A4tech webcam PK -835G, it is a low-cost camera suitable for the prototype because there 
is an adjustment where parameters of the image can be manipulated. The developed prototype in this study is 
shown in Figure 7. 





Figure 6. Prototype bottom view showing the LED Figure 7. The prototype with the height 
strip and webcam adjustment set-up 


The design of the prototype was based on a microscope wherein the height of the camera can be 
adjusted. For most of the testing done in the prototype, the camera module stopper was set-up at 14.5 cm. 
At this height, the camera can take 40 beans at a time. The lighting was optimum at this height without 
producing any shadows. More beans can be captured at higher set-up, but the additional lighting needs to be 
added. The based is from an old electronic drill used for drilling printed circuit board. 


3. RESULTS AND DISCUSSION 
3.1. The list of features that can be extracted by the prototype 

More than 2000 Robusta green coffee beans were used for testing and evaluation of the prototype. 
The green coffee beans are from Indang, Cavite from which morphology, color and texture features were 
extracted. Table 2 shows the list of features that can be extracted by the program. There are a total of 21 features 
that can be extracted by the program. It is composed of 11 morphological features, 6 color features and 4 texture 
features. 


3.2. Comparison of the developed algorithm with the Sobel, Prewitt, and Roberts edge detection algorithm 

The developed edge detection algorithm was loaded in a computer and its performance was compared 
to Sobel, Prewitt and Roberts Edge algorithm as shown in Figure 8. Using visual observation on Figure 8, 
and it can be seen that the developed algorithm is better in detecting edges as compared to Sobel, Prewitt, 
and Roberts algorithm. The images in Figure 8 were compared using the mean squared error (MSE), peak 
signal to noise ratio (PSNR) and execution times. All the images are made same size, because the Python 
program functions only if they are of different sizes. Execution times of the newly developed algorithm 
and the three first order edge detection algorithms were independently timed, measured and compared. 
The result of the evaluation is shown in Table 3. 

The developed algorithm has lower mean square error compared to Sobel, Prewitt and Roberts edge 
detection algorithm which means that the difference between the original image and the image with 
the detected edge is lowest as compared to the other three algorithms. The developed algorithm got the higher 
PSNR, this confirms the visual inspection which means that the peak signal is higher in the developed algorithm 
as compared to noise. Also, execution times in the developed algorithm is faster than all three first degree 
derivative algorithms. 

Better performance of the developed algorithm as compared to other three conventional algorithms 
can be attributed to the adjustments done by normalization to the image. In Sobel, Prewitt and Roberts 
edge detection shown in Figure 8, the noise inside the edge of the coffee beans can be clearly seen, whereas 
in using the newly developed edge detection these noises were eliminated. Also, the execution time of 
the newly developed algorithm is 2.49 seconds which is faster than all three conventional first degree 
derivative algorithms. 


Green coffee beans feature extractor using image processing (Edwin R. Arboleda) 


2032 O 


ISSN: 1693-6930 


Table 2. Features that can be extracted by the program 


Morphology 
Area 
Perimeter 
Length 
Width 


Rectangular Aspect Ratio 


Aspect Ratio 
Major Axis 
Minor Axis 
Roundness 
Ferret Diameter 
Eccentricity 


Color 
Average Red 
Average Blue 
Average Green 
Average Hue 
Average Saturation 


Texture 
Entropy 
Contrast 
Energy 
Homogeneity 


Average Value (Brightness) 


Table 3. The developed algorithm Vs. Sobel Vs. Prewitt Vs. Roberts 


Edge Detector MSE 
This study 68. 55 
Sobel 78. 85 
Prewitt 74. 10 
Roberts 80.95 
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©0000 000 
Q 0000000 
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PSNR 
29.11 
29.16 
29.43 
29.05 










Execution Times (S) 
2.49 
8.04 


7.64 
1.21 
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Figure 8. Original image and the different edges derive using different algorithms; (a) original image, 
(b) newly developed edge detection algorithm, (c) Sobel edge detection algorithm, 
(d) Prewitt edge detection algorithm, (e) Roberts edge detection algorithm 
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4. CONCLUSIONS 
Based on the results of the different tests conducted in this study, the following conclusions are drawn: 
A discriminant for edge detection was developed as defined as 


28 ee l Intensity 
New Edge Detection Discriminant = Total Difference « ————__—__ 
Normalization 


As applied to coffee beans, the optimal values are 10 for intensity; and 25 to 40 for normalization 
and dependent on image. The pixel pair formation that produces optimal peak signal to noise ratio is the sum 
of the difference of top pixel and bottom pixel and difference of left pixel and right pixel. The improvement of 
the developed edge detection algorithm over Sobel, Prewitt, and Roberts are: 2.091%, 1.16%, and 2.47% in 
terms of PSNR, respectively; 13.06%, 7.48%, and 15.31% in terms of MSE, respectively; 69.02%, 67.40%, 
and 65.46% in terms of execution time, respectively. The coffee bean features extracted by 
the prototype are: area, perimeter, length, width, rectangular aspect ratio, aspect ratio, major axis, minor axis, 
roundness, ferret diameter, eccentricity for morphology; average red, average blue, average green, average hue, 
average saturation, average values for color; and entropy, contrast, energy, and homogeneity for texture. 

Overall, it can be claimed that the general objective of this study which is to develop a new edge 
detection approach for green coffee beans feature has been achieved. The main objective has been achieved by 
using the heuristic approach in calculating the right values for the discriminant and finding the best pixel 
formation that produces the optimal PSNR. The newly developed edge detection algorithm was proven better 
than Roberts, Prewitt and Sobel algorithm in terms of PSNR, MSE and execution times. Upon integrating 
the newly developed algorithm in a Raspberry Pi microcomputer, the developed hardware was able to extract 
the morphological, color and texture features of the green coffee beans. 
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